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Abstract 

In this paper, we review the roles of collective modes in nuclear re- 
actions. We emphaize the strong couplings of various collective states 
with the monopole and quadrupole motions. In inelastic excitation, 
these couplings can be seen as an important source of anharmonicity 
in the multiphonon spectrum. In fusion reaction, the breathing and 
quadrupole motions strongly affects the oscillation of protons against 
neutrons. Finally, the modification of the collective properties in- 
duced by a large amplitude dilution might be the origin of the nuclear 
multifragmentation, directly related to the nuclear liquid- gas phase 
transition. In the three cases we derive a coupling matrix element 
which appears to be in good agreement. 

1 INTRODUCTION 

Strongly interacting systems with many degrees of freedom are the prototypes 
of complex systems. As a consequence of this complexity, their dynamics is 
expected to present disorder or chaos. However, in reaction to an external 
stress, such systems appear to often self-organize in simple collective motions. 
Of particular importance is the occurrence of collective vibrations which, in 
general, are surprisingly harmonic despite their chaotic environment. 

This paradox is well illustrated by the atomic nucleus. Indeed, on the 
one hand, following the Bohr ideas, the compound nucleus resonances sign 
the occurrence of quantum chaos already just above the neutron threshold 
p. On the other hand, in the same excitation energy domain, the nucleus 
is known to exhibit a large variety of collective vibrations (called phonons) 
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[2]. The first quantum of oscillation is associated with the giant resonances, 
anomalously large cross sections observed in some nuclear reactions. The first 
one to be discovered was the giant dipole resonance (GDR) [3| interpreted 
as a collective oscillation of the neutrons against the protons Then came 
the giant quadrupole resonance (GQR), an oscillation of the nucleus shape 
between prolate to oblate deformations 5J, and the giant monopole resonance 
(GMR) or breathing mode, an alternation of compressions and decompres- 
sions of the whole nucleus jB]. Since then many other resonances have been 
uncovered OIZI- Twenty years ago, giant resonances have also been observed 
in hot nuclei up to several- MeV temperature |HI. This demonstrates the sur- 
vival of ordered vibrations in very excited system which are known to be 
chaotic. 

The study of this amazing self-organization of the nucleus in collective 
vibrations and its transition from order to chaos is one of the important 
subjects in modern nuclear physics. 

In this article, we will focus on the onset of disorder through the coupling 
between various modes. First we will show that collective vibrations induce 
monopole (GMR) and quadrupole (GQR) oscillations. This means that the 
coupling matrix elements between a quantum of vibration and a state with 
a GMR or a GQR built on top of it are large. We will present results 
from two "orthogonal" approaches, i.e. boson mapping (BM) and time- 
dependent mean-field (TDMF), showing the same effect. Then we will move 
to larger amplitude motion reached in reactions. First we will discuss the 
effect of the large amplitude monopole and quadrupole oscillations induced 
by fusion reaction on the GDR. Then, we will move to even more violent 
reactions for which a rapid expansion of the produced nuclear system have 
been observed and interpreted as the result of fast decompression of the 
matter. We will show how this large-amplitude breathing mode affects all 
the other collective states, which may even become unstable. We will also 
make the bridge between this coupling of collective motions and the liquid- 
gas phase transition. 

Finally we will connect the three studied phenomena with a non linear 
coupling between a vibration and a GMR or GQR built on top of it. 
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2 Theoretical framework 



In quantum mechanics, harmonic oscillations are associated with boson de- 
grees of freedom. From the microscopic point of view, these bosons can be 
understood as being built from fermion pairs, which carry boson quantum 
numbers. However, the number of possible pairs must be large enough to 
insure that the effects of the fermion antisymmetrization do not introduce 
significant deviations from a boson behavior. In particular, the excitations 
of small fermionic systems are not expected to be well described by a boson 
picture, because the Pauli exclusion principle imposes constraints that can- 
not be easily accounted for in a boson representation. From a formal point of 
view the relation between fermion pairs and bosons can be explicitly worked 
out using boson mapping techniques |9J. We will use one of these methods 
in the first study we are presenting. 

Fermionic approaches can also be followed. For example, giant resonances 
are often described using time dependent mean field approaches like the 
Time Dependent Hartree-Fock approximation (TDHF). Indeed, they corre- 
spond to the response of the system to an external (collective) one-body 
field and mean-field approaches are tailored to take care of such excitations. 
Moreover, giant resonances directly affect the time evolution of one-body 
(collective) observables which are well predicted by mean-field approaches. 
The small amplitude reduction of these approaches is equivalent to the ran- 
dom phase approximation (RPA) 9J. Being a linearization of the equation of 
motion, it corresponds to a harmonic picture. However, since the mean-field 
depends upon the actual excitation, TDHF is a non linear theory and hence 
contains couplings between collective modes. For example the quadratic re- 
sponse takes into account the couplings between one and two phonon states 
coming from the 3-particle 1-hole and 1-particle 3-hole residual interaction. 
In fact TDHF is optimized for the prediction of the average value of one 
body observables. Through non-linearities and time dependence, it takes 
into account the effects of the residual interaction as soon as the considered 
phenomenon can be observed in the time evolution of a one body observable. 
This was already the case for the RPA, which through the time dependence 
takes into account the particle-hole residual interaction and goes beyond the 
static mean field which is limited to the hole-hole terms. In this article, we 
will go beyond the RPA treatment either working out the quadratic response 
to TDHF or directly performing full TDHF calculations [TU] . 
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3 Multiphonon anharmonicities due to the 
coupling with GMR and GQR 



Let us start with a direct manifestation of the couphng between collective 
states: the anharmonicity of multiphonon spectra. 

3.1 TDHF picture 




monopole boost quodrupole boost dipole boost 



Figure 1: Evolutions of the monopole, quadrupole and dipole moments (solid 
lines) as a function of time for monopole (a), quadrupole (b) and dipole (c) 
excitations in '^°Ca. 

Time dependent approaches provide an intuitive understanding of collec- 
tive motions [TI]. Let us for example look at the TDHF dynamics for the 
^°Ca nucleus which has been initially perturbed by a collective boost. For 
this simulation we have used the TDHF 3D code developed by P. Bonche 
and coworkers ITU] with the SGH Skyrme force |12j . 

In Fig. m we followed the monopole, quadrupole and dipole response for 
three initial conditions: 

• A monopole boost using Qq = — {rf){t = 0)) as a boost 
generator. Because of the spherical symmetry, a monopole boost can 
only trigger monopole modes. Therefore, we only observe {Qo){t). 

• A quadrupole boost generated by Q2 = J^i^f'^o i^i^Vi) ■ The parity 
conservation forbids any dipole excitation when a quadrupole velocity 
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field is applied to a spherical nucleus. Conversely, breathing modes can 
be triggered by the quadrupole oscillation so that we do follow both 
the quadrupole {Q2){i) and the monopole {Qo)(t) responses. 

• An isovector dipole boost induced by = Z/AJ2n^n — N/AJ2pZp. 
This excitation can be both coupled to the quadrupole and monopole 
oscillations so that we monitor the three moments, {QQ)(t), {Q2){t) and 

{Qnm. 

3.1.1 Linear response and collective states 

In Fig. n we observe that the collective boosts induce oscillations of the 
associated moments as expected from the RPA picture. They are only slightly 
damped in the GQR and GDR cases (fig. 1-b and 1-c respectively) while 
in the GMR case (fig. 1-a) beatings, characteristic of a Landau damping 
[TT] . are observed. This means that the dipole and quadrupole strengths 
are mostly concentrated in a single resonance while the monopole one is 
fragmented. 




monopole boost quadrupole boost dipole boost 



Figure 2: Monopole, quadrupole and dipole spectra obtained through the 
Fourier transform of the time dependent response for a monopole, quadrupole 
and dipole excitation in the "^^Ca. 

Looking at the amplitude of the first oscillation {Qv)max as a function 
of boost strength iky) confirms the linearity of this response ^H]. To get a 
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deeper insight into the response, we study the Fourier transform F^{uj) of 
{Qv)it)/K which is nothing but the RPA response when the velocity field 
is small enough to be in the linear regime. We see in Fig. |21that the dipole 
and quadrupole modes are concentrated in a unique peak. The monopole is 
fragmented but the various peaks are in the same energy region so that they 
can be approximated by a single mode with a large Landau width. 

3.1.2 Quadratic response 

In Fig. [T[ one can see that large amplitude dipole (fig. 1-c) and quadrupole 
(fig. 1-b) motion induces variations of the nucleus radius {Qo). Looking at 
the matter distribution one can see that the central density po is affected by 
the dipole or quadrupole motion. Since the central density can be modified 
only by monopole states this imposes that the large amplitude motion gets 
coupled with such breathing modes. In the same way a large amplitude 
dipole oscillation induces a quadrupole deformation of the nuclear potential 
and so gets coupled with the GQR. These observations lead to the conclusion 
that we are in the presence of a non-linear excitation of a giant resonance 
(GMR or GQR, generically called /i in the following) on top of the collective 
motion (GQR or GDR, generically called u in the following) initially excited 
through the collective boost Qi,. As expected from the quadratic response 
theory [T3], the amplitude of the first oscillation of {Qfj,)(t) is quadratic in the 
excitation velocity and its time dependence corresponds to (cos(co'^t) — 1), 
i.e. it is initially in phase quadrature with {Qi,){t) but it is oscillating with 
the frequency of the mode /i and not the one of the initially excited collective 
state z/. 

This interpretation is confirmed by the Fourier transform of (Q^)(t) as- 
sociated with the excitation of Qu (Fig. Let us first start with the 
quadrupole strength non-linearly excited by a dipole boost. The observed 
peak is identical to the GQR-response. This is a clear indication that the 
observed state is indeed a GQR built on top of the GDR. It should be noticed 
that this frequency is different from the one of the underlying dipole motion. 
The monopole case is more complex because of the presence of a strong Lan- 
dau spreading and it seems that the strengths of the various monopole states 
depend upon the considered boost. This indicates that the coupling leading 
to the excitation of an additional monopole state depends upon the collective 
mode initially excited. 
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IjJy 




(z/[V[z/0) 


{v\V\v2) 




(MeV) 




(MeV) 


(MeV) 


0)40Ca 


22.9 


11.6 







2)40Cq 


18.6 


21.4 


-4.28 




\D)AOCa 


17.2 


3.47 


-4.58 


-3.92 


0)208p;, 


15.7 


57.1 







[2)208 Pft 


11.1 


99.0 


-2.17 




[Z))208 PI) 


13.0 


8.94 


-2.40 


-0.70 



Table 1: Energies, transition probabilities and coupling coefficients of the 
GMR, GDR and GQR in the ^°Ca and ^ospb. is expressed in fm^ for the 
GMR and GQR and in fm for the GDR. 



3.1.3 Quadratic response and Couplings between states 

Assuming for each multipolarity a unique state [/i) non linearly excited one 
can use the quadratic response theory 13j to extract the residual interaction 
matrix element between \v) and [i//i) from the amplitudes of the induced 
oscillations {Q^,)raax using 

^ 2klqiq, 

where is the transition matrix element between the ground state [— ) and 
the collective state [77), = {—\Q-q\v)y which can be derived from the linear 
response with t] = u since 

max /„\ 

The results for the ^°Ca and ^°*Pb are presented in table 1. The Ui, are 
computed from the time to reach the ffist maximum of {Qp){t). The relative 
sign of and q^ is given by the early evolution of the moments {Q^) (see 
Eq. They appear to be all negative. The couplings are large of the 
order of few MeV. As we will see those findings are in qualitative agreement 
with the results of ref. ^3] which is using a completely different approach. 
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T 


Eharm 


EWSR 


{y\V\uQi) 


{y\V\uQ2) 


{u\V\u2) 






(MeV) 


(%) 


(MeV) 


(MeV) 


(MeV) 


|Ol)40Ca 


0+ 


18.25 


30 


-2.13 


-2.36 




02)4()Ca 


0+ 


22.47 


54 


-2.03 


-3.96 




ll)4()Ca 


1- 1 


17.78 


56 


-1.38 


-2.12 


-1.25 


ll)4()Ca 


1- 1 


22.03 


10 


-1.48 


-2.16 


+0.73 


2)40Ca 


2+ 


16.91 


85 


-1.36 


-2.49 


-0.36 



Table 2: RPA one-phonon basis for the nucleus ^°Ca. For each state, spin, 
parity, isospin, energy and percentage of the EWSR are reported in addition 
to the coupling coefficients. 





JTT rj. 


E 


EWSR 


{y\v\vQi) 




{u\V\u2) 






(MeV) 


(%) 


(MeV) 


(MeV) 


(MeV) 


0l)2OSpb 


0+ 


13.61 


61 


-1.87 


-0.92 




02 ) 208pfe 


0+ 


15.02 


28 


-1.32 


-1.16 




ll)208p6 


1- 1 


12.43 


63 


-0.79 


-0.59 


-0.68 


12)208 pf, 


1- 1 


16.66 


17 


0.00 


0.00 


-0.64 


2) 208p6 


2+ 


11.60 


76 


-0.64 


-0.48 


-0.74 



Table 3: Same as table for the nucleus ^'^^Pb. 
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3.2 Comparison with boson mapping calculations 



In ref. ^3] a completely different approach is used to infer the same matrix 
elements: the fermionic Hamiltonian is first mapped into a bosonic one mak- 
ing a connection between any particle-hole excitation and a boson. Because 
of the fermionic anticommutation relations (Pauli principle) a particle-hole 
excitation operator is mapped into an infinite series of boson operators. As a 
consequence, even if the fermionic Hamiltonian is containing only two body 
interaction, the boson Hamiltonian is a infinite series with many boson inter- 
actions not conserving the boson number. To be manageable this series has to 
be truncated and, in the application presented in only the terms contain- 
ing up to four boson creation or anihilation operators have been conserved. 
Then a RPA transformation is applied, introducing collective phonons which 
optimizes a harmonic picture to the system properties J3] . Using those col- 
lective degrees of freedom, the Hamiltonian then contains an harmonic part 
plus various interactions. An important one is the coupling between one- 
and two phonon states. 

For practical reasons only the most collective phonons have been selected 
in ref. ^3] some of them are presented in table 2 and 3 for ^°Ca and ^°^Pb. 
To illustrate their degree of collectivity the part of the energy weighted sum 
rule (EWSR) they are exhausting is given together with their energy. Table 
2 and 3 also give the coupling matrix elements between different collective 
states and one of the two GMRs or the GQR built on top of them. 

From the quantitative point of view, the non-linear coupling extracted 
from TDHF appears to be 50% larger than the one reported in Table 2 
and 3 ^ll- This is a reasonable agreement since TDHF re-sums all the 
individual couplings as shown in |13j. Summing the contributions of the 
different collective states considered in ref. [H| reduces the difference between 
the reported values. However, the phonon basis studied in ref. ^3] being 
incomplete , it is expected that the TDHF results remains higher. It should 
be also noticed that some differences can remain due to the approximations 
involved in the different approaches. 

3.3 Consequences on the multiphonon spectrum 

The couplings between phonon states can be used to derive the energies of the 
various phonon states. It appears that the very large matrix elements exciting 
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a GMR or a GQR on top of any state are the main source of anharmonicites. 
As far as the two-phonon states are concerned the couphng with the three- 
phonon states coming from such a non-hnear couphng induces an energy shift 
of about 2 MeV for Ca and 1/2 MeV for Pb. This might be the explanation 
of the phonon anharmonicity puzzle discussed by many authors jl6[ [T7| ITH] 
and summarized in 

4 Pre-equilibrium GDR in fusion reactions 

It has been recently proposed that GR can play an essential role in fusion 
reactions. In particular, fusion reactions with N/Z asymmetric nuclei may 
lead to the excitation of a GDR because of the presence of a net dipole mo- 
ment in the entrance channel ^Hl- Experimental indications of the possible 
existence of such new phenomenon have been reported in refs. fK^ 1^ |22] 
for fusion reactions and [23 121| for deep inelastic collisions. 

4.1 Collective dynamics in fusion. 

As an example ,j25j, we have computed the central collision of Mg 
at 1 MeV per nucleon with the TDHF approach. The system rapidly fuses 
producing an excited ^°Ca nucleus. It presents a strong quadrupole oscillation 
around a slowly damped deformation. Since has a N over Z ratio different 
from this of ^°Mg (respectively 1.5 and 0.67), the center of mass of the protons 
is initially different from the neutrons one. As the time goes on this dipole 
moment Qd (i.e. the distance between the neutron and proton center of 
mass) oscillates (see Fig. 3). To study the induced motion one can plot the 
dipole moment as a function of the velocity of protons against neutrons 
which can be considered as its conjugated moment, P^. We observe a spiral 
in this collective phase space {Q^,PJ , i.e. oscillations in phase quadrature 
of the conjugated dipole variables. This is a clear signal of the presence of a 
damped collective vibration. 

The period of the observed oscillations is around 150 fm/c while for the 
^°Ca nucleus in its ground state it is almost half this value. This large dif- 
ference can be explained by the deformation of the fused system. Indeed, 
in the TDHF simulations the compound nucleus only slowly relaxes its ini- 
tial prolate elongation along the axis of the collision. The averaged value 
of the observed quadrupole deformation parameter is around e = 0.23. For 
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Figure 3: Time evolution of the dipole vibration. Dipole moment Qd and its 
conjugate are plotted in the phase space (a) and Pd is plotted in function 
of time (b). 

symmetry reasons, the dipole oscillation occurs only along the deformation 
axis of the nucleus formed by fusion in head on head reactions. Therefore, 
a lower mean energy is expected for this longitudinal collective motion ac- 
cording to the following relation Eq£,r^ = Egdr (1 ~ The energy of the 
GDR along the elongation axis Egdrz fuUfils this relation with £ fti 0.26 in 
a good agreement with the previous calculation of e. 

4.2 Dynamical coupling of the GDR with the nucleus 
deformation 

To get a deeper insight in the dipole oscillation observed in fusion reactions 
we have analyzed the time evolution of the period. From each point on the 
collective trajectory in the collective phase space {Q^,P^) this quantity can 
be inferred from the time needed to reach the opposite side of the observed 
spiral (see Fig. 3). The resulting evolution is plotted in Fig. 4 — a. This 
period presents oscillations too. These variations of the GDR period are 
almost in phase with the observed oscillations of the monopole moment Qo 
and the quadrupole moment Q2 presented in Fig. 4 — 6. This points to a 
possible coupling between the dipole mode and another mode of vibration. 
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The evolutions of the monopole and quadrupole moments are very similar. 
In particular, they present the same oscillation's period around 166 fm/c. 
Therefore, they originates from the same phenomenon, the vibration of the 
density around a prolate shape. This oscillation modifies the properties of 
the dipole mode in a time dependent way. 



^240 




Monopole 



100 200 300 400 500 600 700 300 

Time (fm/c) 



Figure 4: (a) Time behaviour of the dipole period (solid line) and its mod- 

elisation with the Mathieu's equation (dashed line), (b) Time evolution of 
the monopole (dashed line) and quadrupole moments (solid line). 

In order to investigate if this is the origin of the observed behavior we can 
model the induced mode coupling. Let us consider a harmonic oscillator of 
average frequency c^o and let us replace the influence of the density oscillation 
characterized by the frequency uj by assuming that its spring constant varies 
in time at a frequency uj. In such a model, the equation of the motion 
becomes the Mathieu's equation 

^ + [1 + 5cos(cjt)]x = (3) 

where c^o is the pulsation without coupling and oj the pulsation of the den- 
sity's oscillation while 5 corresponds to the magnitude of the induced fre- 
quency fluctuations. We have computed the numeric solution of Eq. 01 with 
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the typical external frequency of the monopole and quadrupole oscillations. 
The bare frequency Uq and the coupling strength S have been tuned in order 
to get the same oscillation frequency and typical amplitude. Indeed, because 
the presence of the oscillating term the observed frequency is different from 
the bare one. The solution (see Fig. 4-a) well reproduces our observations 
with 5 — —0.32 and cuq changed by a factor 1.2 from the observed value 
EgdRz- From this analysis it appears that the observed dipole motion cor- 
responds to a giant vibration along the main axis of a fluctuating prolate 
shape. 

4.3 Couplings between modes 

It is interesting to relate 5 to a coupling matrix element between collective 
states. Indeed, the Mathieu's equation can be seen as the equation of mo- 
tion coming from an Hamiltonian containing a coupling term between the 
dipole mode and the collective deformation (GMR or GQR called here after 
ii)V= I^TT^-r — QnQ^ ■ Thls leads to a coupling matrix element between the 

GDR and the state ii built on top of it which reads V = S^-n=r^ — B\b\B,. 
In the studied case for both the monopole and quadrupole, using the ground 
state matrix elements g^, the amplitude of the oscillation ^'^^J'""'^ is about 5. 
Using a;o — 13 MeV we get a coupling between the dipole and the monopole 
or quadrupole v — S^j^^ — — —0.5 MeV. This value is qualitatively in 
agreement with the previous observation of a strong coupling. From a quan- 
titative point of view, it is lower than the corresponding one (see table 1). 
However, it is obtained for a hot and deformed system with a large ampli- 
tude monopole and quadrupole motion. Moreover the value derived here 
only correspond to the coupling with a unique mode and so should rather be 
compared with table 2. 

5 GR in diluted nuclei 

In violent heavy ion collisions, nuclear matter is excited and compressed. 
Then the formed nuclear system expands under the resulting thermal and 
mechanical pressure. Matter may also be quenched in the coexistence region 
of the nuclear liquid-gas phase diagram and the observed abundant frag- 
ment formation may take place through a rapid amplification of spinodal 
instabilities. New experimental results pleading in favor of such a spinodal 
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decomposition have recently been reported [^Hl 123 • From the theoretical 
point of view, the spinodal instabilities in finite systems are unstable col- 
lective motions. They have been mainly studied within semi-classical or 
hydrodynamical framework [2HI 123 EOl El E21 • However, since the relevant 
temperatures are comparable to the shell spacing and the wave numbers of 
the unstable modes are of the order of Fermi momentum, quantum effects 
are expected to be important. Quantum approaches linking the spinodal 
instabilities with the giant resonances can be found in ES] • 

5.1 RPA in diluted systems 

To investigate instabilities encountered during the evolution of an expanding 
system, one should study the dynamics of the small deviations Sp{t) around 
the TDHF trajectory p{t) ISE]- It is more convenient to carry out such an 
investigation in the " co-moving frame" and if we consider the early evolution 
of instabilities in the vicinity of an initial state the problem reduces to an 
RPA like equation . Then, small density fluctuations are characterized by 
the RPA modes pu and the associated frequencies Uu- When the frequency 
of a mode drops to zero and then becomes imaginary, the system enters an 
instability region. 

In order to perform an extensive study of instabilities we may parametrize 
the possible densities po either by a static Hartree-Fock (HF) calculations 
constrained by a set of collective operators ^7\, or using a direct param- 
eterization of the density matrix. In the following, we follow the second 
approach by introducing a self-similar scaling of the HF density as suggested 
by dynamical simulations. 

First, we solve the HF equation for the ground state [hnr^PHF] = 0, 
leading to the single-particle wave functions \ipi) and the associated ener- 
gies Si- Then we introduce the density matrix at a flnite temperature T as 
Phf [T] = 1/ -I- exp [{hnF-^F f'^))' where ep [T] is the corresponding 
Fermi level that is tuned in order to get the correct particle number. Then, 
we perform a scaling transformation, R [a] , which inflates the wave functions 
in radial direction by a factor a according to {r\R [a] \ip) = a~^^^ {r / a\(f) . We 
then deflne the density matrix for a hot and diluted system by po [(^,T] = 
R [a] Phf [(^"^T] R^ [a]. The eigenstates of the constrained Hamiltonian are 
given by \i) = R [a] \ipi), and the corresponding energies and occupation num- 
bers are = Ei/a'^ and = 1/ f 1 + exp ((si-ep [a^T]j /a^Tj j, respectively 
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(see Ref. jSH] for more details). 




Figure 5: Contour plots of the isoscalar strength functions associated with 
the multipolarity L = 2 — 5 as a function of the dilution parameter a and 
the collective energy of the mode = cUujy (c = 1 for stable modes, —i for 
unstable modes) for ^''Ca (left) and ^^°Sn (right). 

We performed the HF calculations in the coordinate representation using 
the Skyrme force SLy4 [2H1- The particle and hole states are obtained by 
diagonalizing the HF Hamiltonian in a large harmonic oscillator representa- 
tion [39 , which includes 12 major shells for Ca isotopes and 15 for Sn. We 
apply the scaling and heating procedures described above to the density ma- 
trix, and calculate the residual interaction in a self consistent manner. We 
solve the RPA by a direct diagonalization using a discrete two quasi-particle 
excitation representation jlUj . 



5.2 Dilution-dependent GR frequencies 

The top part of Fig. |31 shows calculations performed for ^°Ca. Top panels 
shows contour plots of the isoscalar strength function associated with the 
isoscalar operator A^^j^ = J2iLi''"i'^LM, with multipolarity L = 2 — 5, as a 
function of the dilution parameter a. We observe that, in the stable domain, 
the energy associated with the dominant isoscalar strength decreases as di- 
lution becomes larger, and at a critical dilution it drops to zero. At larger 
dilution, the system becomes unstable, and for each multipolarity, one or 
two unstable modes appear. This is illustrated in the bottom panel, where 
the "imaginary energy" of the mode E^, = —ihuOy is plotted as a function 
of the dilution. Looking at the RPA solution in the coordinate space, it 
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appears that the collective motions transform into volume modes when they 
become unstable. We observe in fact a quite complex structure of the un- 
stable modes: Volume and surface instabilities are generally coupled, as well 
as isoscalar and isovector excitations, since protons and neutrons move in 
oposite way jH^j . 

5.3 RPA-instabilities "phase diagram" 

We, also, carry out calculations at finite temperature and determine the 
dilutions at which different unstable modes begin to appear. This allows 
us to specify the border of the instability region in the density-temperature 
plane for different unstable modes. Fig. 4 shows phase diagrams for octupole 
instabilities in ^^''Sn. Here, for simplicity, we define the density a.s p = po/a^. 



— ' — I — ' — I — ' — 
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Figure 6: Border of the instability region (full fine) associated with L = 3, 
for ^^''Sn. The dashed line connects the points having the instability growth 
time r = 100 fm/c. The dots are associated with r = 50 fm/c. 

The full line indicates the border of the instability region. The dashed line 
connects points that are associated with the instability growth time r = 
100 fm/c, and the dots correspond to situations with a shorter growth time 
r = 50 fm/c. We observe that in finite nuclei the instability region is quite 
reduced as compared to that of nuclear matter. The limiting temperature 
for instability to occur is around 6 MeV for Sn and 4.5 MeV for Ca (see |35j ) 
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while it is about 16 MeV in symmetric nuclear matter. As a result, heavier 
systems have larger instability region than the lighter ones. 



5.4 Link with the coupUng between modes 

The observed dilution dependence of the GR's energies can be interpreted 
in terms of a coupling between the studied GR and the GMR. Indeed, it 
corresponds to the Hamiltonian H = uj [a) B\Bj where uj [a) is the GR fre- 
quency for a dilution a. Using a Taylor expansion of ui (a) around a = 1 
and introducing e — a — 1, cu (a) = a;o -l- edaUjlg, and the Hamiltonian 
becomes H — uqB\B^ -\- e da,oj\Q B\B ^. The dilution factor a can be re- 
lated to the collective observables using (r^) = a^(r^)o — {r'^)o + 2£(r^)o 
and — (r^ — (r^)o) = % (^Bj^ + B^^j so that the Hamiltonian reads 

H — ujqB\B^ -\- l^y^ daOj\QB^^B\B^ -\- h.c. . This leads to a couphng ma- 
trix element between the GR and a GMR (//) built on top of it which reads 
y — "^py^ ^a^loBj^Bl^B^. Looking at the GQR decrease with the dilution, 

we get daUJ ~ —40 McV for ^''Ca. Using the transition matrix clement for 
the GMR extracted in table 1, ~ 12 fm^, and (r^)o/A/47r ~ 125 fm^ we 



get Vf^ = 2(t-2)^7v^ 9aUj\Q ~ —2 McV. This value is qualitatively in agreement 
with the previous observation of a strong coupling. In particular we under- 
stand the sign of the interaction since the energy of the modes is reduced in 
a diluted system < 0). From a quantitative point of view, it is a factor 
2 lower than the one reproduced in table 1. However one should remember 
that we compute here an average matrix clement while in table 1 it is more 
an integrated one. In fact this result should be compared to the one in table 
2 and indeed they are in excellent agreement. 



6 Conclusion 

In conclusion, we have presented a study of 3 different situations: the excita- 
tion of small amplitude motions and the coupling between phonon states, the 
excitation of pre- equilibrium GDR in fusion reaction and the modification of 
the GDR properties due to the dynamics of the monopole and quadrupole 
deformation and finally the early development of spinodal instabilities. The 3 
phenomena show that collective modes are strongly coupled to the monopole 
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(and quadrupole) degrees of freedom. We have presented original methods to 
extract the couphng matrix element between a collective mode and a GMR 
(or GQR) built on top of it. The 3 studied phenomena are in qualitative 
agreement pointing to a negative coupling of the order of the MeV. The 
quantitative difference we observe which might come the differences between 
the 3 studied cases is now under study. 
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